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ABSTRACT 

Using a samples of 61506 spiral galaxies selected from the SDSS DR2, we examine the luminosity 
function (LF) of spiral galaxies with different inclination angles. We find that the characteristic 
luminosity of the LF, L*, decreases with increasing inclination, while the faint-end slope, a, depends 
only weakly on it. The inclination-dependence of the LF is consistent with that expected from a 
simple model where the optical depth is proportional to the cosine of the inclination angle, and we use 
a likelihood method to recover both the coefficient in front of the cosine, 7, and the LF for galaxies 
viewed face-on. The value of 7 is quite independent of galaxy luminosity in a given band, and the 
values of 7 obtained in this way for the 5 SDSS bands give an extinction curve which is a power 
law of wavelength (r oc A""), with a power index n = 0.96 ± 0.04. Using the dust extinction for 
galaxies obtained by Kauffmann et al. (2003), we derive an 'extinction-corrected' luminosity function 
for spiral galaxies. Dust extinction makes M* dimmer by about 0.5 magnitudes in the z-band, and 
about 1.2 magnitudes in the u- band. Since our analysis is based on a sample where selection effects 
are well under control, the dimming of edge-on galaxies relative to face-on galaxies is best explained 
by assuming that galaxy disks are optically thick in dust absorptions. 
Subject headings: galaxies: spiral — galaxies: luminosity function — ISM: dust, extinction 



1. INTRODUCTION 

Dust plays an important role in the observed properties 
of galaxies. It not only causes galaxies to dim but also 
causes them to redden. In order to obtain the intrinsic 
properties of galaxies from observation, it is necessary to 
understand dust extinction. Dust extinction is believed 
to be more important in late-type spiral galaxies than 
in early- type galaxies, since spiral galaxies are richer in 
cold gas and have on-going star formation (e.g. Calzetti 
2001 for a review and references therein). Dust opacity 
in spiral disks can be probed by studying the average 
photometric properties of spiral galaxies as a function 
of disk inclination angle, as first suggested by Holmberg 
(1958, 1975). Theoretically, it can be shown that, for 
completely opaque disks, the surface brightness to be ob- 
served is almost independent of inclination angle, while 
the luminosity dims as the inclination angle changes from 
face-on to edge-on. On the other hand, for completely 
transparent disks, the luminosity is independent of incli- 
nation while the surface brightness brightens as the disk 
changes from face-on to edge-on. However, the power of 
the Holmberg test is limited by the finite sizes of galaxy 
samples, because it is in general difficult to accurately 
estimate the inclination angles for individual galaxies. 
The axis ratios of galaxy images are usually used to es- 
timate the inclination angles. For an infinitesimally thin 
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and round disk, the inclination angle 6 can be obtained 
directly from the axis ratio b/a through cos9 = b/a. Un- 
fortunately, real disks are neither completely round nor 
infinitesimally thin, and the observed axis ratio depends 
not only on the inclination angle but also on the elliptic- 
ity and thickness of the galaxy. Unless such dependence 
is fully taken into account, the results from the Holm- 
berg test cannot be interpreted straightforwardly. In ad- 
dition, the test based on the average photometric prop- 
erties of galaxies as functions of inclination angle may be 
further complicated by the incompleteness and selection 
bias of galaxy sample. For example, if disk galaxies are 
optically thin, a galaxy sample may be biased for low- 
luminosity galaxies with high inclinations because of the 
the enhanced surface brightness. With such selection 
bias, face-on galaxies will on average be brighter than 
edge-on galaxies, which may be falsely interpreted as the 
disks being optically thick. Because of these limitations, 
the conclusion about the opacity of the disk galaxies are 
still controversial (e.g. Burstein Haynes & Faber 1991; 
Byun 1993; Giovanelh et al. 1994, 1995; Xilouris et al. 
1999; Masters, Giovanelli & Haynes 2003; Holwerda et 
al. 2005). 

With a large and well-defined sample, the above- 
mentioned difficulties in the Holmberg test can in princi- 
ple be overcome. As shown in Ryden (2004), a large sam- 
ple allows one to estimate the thickness and ellipticity ro- 
bustly from the distribution of axis ratios in a statistical 
way. Moreover, with a large, well-defined sample, one can 
study the conditional distribution functions of photomet- 
ric properties for given inclination angles, instead of just 
considering the mean photometric properties of galaxies. 
Thus, the selection bias and sample incompleteness can 
be taken into account strictly in a statistical sense. With 
the advent of the Sloan Digital Sky Survey (SDSS, York 
et al. 2000), it is now possible to make such analysis. 
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In this paper, we use the SDSS galaxy sample to study 
dust extinctions in spiral galaxies. Our analysis is based 
on the luminosity function (hereafter LF) of galaxies as 
a function of inclination angles. We statistically correct 
the axis ratios of galaxy images to obtain their inclination 
angles based on the method of Ryden (2004). We com- 
pare the change of LF with inclination angle in different 
wavebands to constrain the shape of the dust extinction 
curve. The paper is organized as follows. The selection 
of sample is described in § 2.1 and quantities to specify 
the inclinations of galaxies are discussed in § 2.2. In § 3, 
after a brief description of the LF estimators adopted in 
this paper (§ 3.1), we present the results of the LFs in 
all the 5 SDSS bands for galaxies of different inclination 
angles (§ 3.2). In § 4, the inclination-dependence of the 
LF is modelled in terms of dust extinction. The possi- 
ble luminosity dependence of the dust extinction and the 
dust-corrected LF of spiral galaxies are discussed in § 5. 
Finally, our results are summarized in § 6. 

2. OBSERVATIONAL DATA 

2.1. Galaxy samples 

The galaxy samples we used were selected from 
the New York University Value-Added Galaxy Catalog 
(NYU-VAGC, Blanton et al., 2005) of the SDSS second 
data release (DR2, Abazajian et al. 2004). The red- 
shift catalog of the DR2 covers 2,627 deg^ of the celestial 
sphere and photometric data in the 5 SDSS wave bands, 
u, g, r, i and z, are available for each of the galaxies 
(Abazajian et al. 2004) directly from the SDSS pipeline. 
The NYU-VAGC includes additional information for ex- 
tragalactic targets, such as if-correction, spectroscopic 
target completeness, etc.. From the NYU-VAGC, we se- 
lect spiral galaxies with fracdeVr < 0.5, where the pho- 
tometric parameter fracdeV is a point spread function 
(PSF) corrected indicator of galaxy morphology. In the 
SDSS pipeline, each galaxy was fitted by an exponential 
profile and a de Vaucouleurs' profile. The best linear 
combination of these two profiles was used to represent 
the profile of the galaxy, and fracdeV is the fraction 
of luminosity contributed by the de Vaucouleurs' profile. 
Bernardi et al. (2005) used fracdeVr > 0.8 to select 
early-type galaxies. We use fracdeVr < 0.5 to ensure 
that the galaxies we select are dominated by the expo- 
nential component. 

To construct samples to study the LFs of galaxies, both 
redshift and flux limits are applied. Only galaxies with 
redshift in the range 0.01 < Zrod < 0.22 are selected. 
The upper limit is imposed to minimize the uncertainty 
in the iiT-correction and possible redshift evolutions of 
galaxies; the lower limit is employed to avoid the nearest 
galaxies with large uncertainties of distance caused by 
peculiar velocities. A flux limit in r band, rur = 17.60, is 
chosen to ensure the completeness in spectroscopy. For 
other wave bands, the flux limits are chosen to be suf- 
ficiently high, so that almost all of the galaxies (99%) 
brighter than the limit in the band in consideration have 
r-band flux brighter than 17.60 mag. The flux limits of 
all 5 bands we use and the numbers of galaxies in the 
corresponding samples are listed in Table 1. Note that 
we use the Petrosian (1976) magnitude to refer the lumi- 
nosities of galaxies. This definition of magnitude has the 
advantage that it is a uniform measurement of flux quite 
independent of the distance of the target. Furthermore, 



TABLE 1 

The Sample of galaxies in Each SDSS Band 



Band Flux Limit Redshift Limit Af^-^l A^gpirai 

u 14.50 < m < 18.60 0.01 <z < 0.22 42033 24106 

g 14.50 <m< 17.90 0.01 <z< 0.22 93221 43252 

r 14.50 <m< 17.60 0.01 <z< 0.22 162279 61506 

i 14.50 <m< 17.15 0.01 <z< 0.22 148934 52336 

z 14.50 <m< 16.85 0.01 <z< 0.22 136151 43653 

Afg^l is the number of galaxies selected within the flux and red- 
shift limits. Afgpijj^i is the number of spiral galaxies selected with 
fracdeVr < 0.5. 



for an exponential profile the Petrosian magnitude con- 
tains almost all of the flux of the source (e.g. Stoughton 
et al. 2002). 

For each band, the absolute magnitude M is calculated 
using 

Mx = mx-51og{dL/10pc)-Kx{z) (1) 

where A is the wave band in consideration, nix is the ap- 
parent magnitude corrected for Galactic extinction based 
on the E{B — V) maps of Schlegel, Finkbeiner & Davis 
(1998), dL is the luminosity distance, and Kx{z) is the 
/C-correction value taken from NYU-VAGC. Through- 
out the paper, we use the standard A— cosmology with 
= 0.3, = 0.7 and Hq = 100/i-ikms-iMpc. 

2.2. Inclination parameters 

The apparent axis ratio of the image of a galaxy, which 
is provided by the SDSS pipeline, is a direct measurement 
of the inclination of spiral galaxies. If disks are thin and 
round, the minor-to-major axis ratio b/a is related to the 
inclination angle 6 (defined to be the angle between the 
line-of-sight and the axis of angular momentum of the 
disk) by b/a = cos9. In this paper, we use the r band 
axis ratio ab-cxpr, taken from the best fit of the images 
of galaxies with an exponential profile convolved with 
the PSF (see Stoughton et al. 2002). Since the sample 
galaxies we selected are expected to be dominated by the 
exponential component [fracdeVr < 0.5), the parameter 
ab-cxpr (hereafter referred to as b/a) should be a reason- 
able representation of the apparent axis ratio. 

Figure 1 shows the distribution of b/a for the sample 
galaxies (solid line). If disks are completely round and 
infinitesimally thin, and if their rotation axes are ran- 
domly orientated in space, the distribution of b/a should 
be flat. For our sample, the distribution is almost flat 
for 0.3 < b/a < 0.9, but disks with axis ratios close to 
1 or close to are less common. The lack of very round 
images {b/a ^ 1) is caused by the intrinsic ellipticity, e, 
of spiral galaxies, as demonstrated in detail by Ryden 
(2004). The lack of images with small b/a is caused by 
two factors. The first is the intrinsic thickness of galaxy 
disks: the typical scale height of a disk is about 10% of 
the disk scale length, as shown in Giovanelli et al. (1994). 
The other is the existence of a central bulge which, be- 
ing sphca'oidal, can thicken the image of a spiral galaxy, 
especially an edge-on spiral. In what follows, we do not 
distinguish these two different causes of disk thickening. 
Instead we use a single 'effective thickness' parameter, 
v, defined as the ratio of scale height to scale length, to 
represent the thickness of a spiral galaxy. 
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Fig. 1. — The distributions of the apparent axis ratio hja (solid 
line) and of the simulated inclinations cos Q (dashed line) . 



Using a sample of SDSS spiral galaxies that have rel- 
ative large apparent size, iig > S-^Tpsf , where i?e is 
the effective radius of a exponential profile and Tpgf 
is the adaptive second moment of the PSF image (so 
that \jTpaf may be used as a characteristic scale of 
the PSF), Ryden (2004) found that both v and loge 
have roughly Gaussian distributions. For example, in 
the r-band, Ryden (2004) found that the median (i/q) 
and dispersion (a^) of the i^-distribution are 0.216 and 
0.056, respectively, while those for the log e-distribution 
are logeo = —1.83 and cTioge = 0.93. However, with a 
typical PSF of 1.5 arcsec, the majority of SDSS galax- 
ies have Re < 5^7^ (e.g. Shen et al. 2003). More- 
over, the apparent thickness is expected to vary with 
the apparent size (Re), with smaller images appearing 
rounder, even though the images are de-convolved with 
the PSF. Galaxies with small apparent sizes include not 
only intrinsically small galaxies at small distances, but 
also intrinsically luminous galaxies at large distances. 
Since the average bulge-to-disk ratio of spiral galaxies 
increases with galaxy luminosity, the bulge components 
may become more dominant in a sample that contains 
more distant galaxies with small apparent sizes. Further- 
more the bulge components are in general easier to ob- 
serve than the disk components, because of their higher 
surface brightness. Both factors will cause the distant, 
apparently small spirals to be rounder. For this reason, 
we divide our sample into subsamples according the ap- 
parent sizes (Re) of galaxies, and study the thickness 
parameter separately for each of the subsamples. 

We assume that both the thickness parameter f and el- 
lipticity log e follow Gaussian distributions and estimate 
the thickness parameters vq and a„ for spirals with differ- 
ent apparent sizes using the following Monte-carlo proce- 
dure. For a given galaxy with given shape parameters {ly 
and e) and viewing angles {0, (p), the apparent axis ratio 
b/a can be calculated using 



b/a = 



A + C-^{A-Cf + B 
A + C + ^/{A-C)^ + B 



1/2 



where 



^ = [1 - e(2 - e) sin^ <fi\ cos^ 6 + u'^ sin^ 6, 
B = 4e^(2 - ef cos^ sin^ <p cos^ ip, 



(2) 

(3) 
(4) 



log Rg(arcsec) 

Fig. 2. — The best-fit values of the effective thickness parameters 
fc'o and (7u (shown as error bars) for spiral galaxies with different 
apparent size log Re in the r-band. Dashed line is the distribution 
of log Re for sample galaxies. 



C = 1 - e(2 - e) cos^ if 



(5) 



(see Binney 1985). Model galaxies are assumed to have 
random distribution in the view angle, and to have ran- 
dom distributions in and log e according to their Gaus- 
sian distribution functions. To reduce model parameters, 
we assume that loge and cioge have the values given by 
Ryden (2004). The best values of F and CTj/ are obtained 
by matching the predicted b/a distribution with the ob- 
served b/a distribution (for each i?e-subsample), with the 
use of the least-square criterion. The fitting results for 
fo and a,y are plotted in Figure 2. As expected, disk 
galaxies with large apparent sizes have smaller 'effective 
thickness'. At apparent size i?,. > 8 arcsec, i^o reaches a 
constant value ^ 0.14. This value may be considered as 
the upper limit of the effective thickness of the disks, or 
be explained as the intrinsic thickness of pure disk should 
be less then 0.14. Since the typical seeing condition in 
SDSS has Tpst = 4Pixcl^, the galaxies selected by Ryden 
(2004) have Rg > 5-^7^ and so are dominated by galax- 
ies with Re > 4: arcsec (see Figure 2). Our results for 
galaxies with i?e ^ 4 arcsec are I'o ± ~ 0.195 ± 0.035, 
close to the results of Ryden 0.216 ± 0.056. 

The above procedure also provides a statistical way 
to relate the real inclination of a galaxy (represented by 
cos^) to the apparent axis ratio (b/a). To do this, each 
galaxy is assigned a possible viewing angle using a Monte 
Carlo method. For a galaxy with given (observed) b/a 
and Re , we first randomly select a set of shape parame- 
ters (i/ and e) from their Gaussian distributions obtained 
above, and then calculate the apparent axis ratio using 
equation (2). The value of that best reproduces the 
observed b/a is taken to be the viewing angle of this 
galaxy. Since our samples are sufficiently large, the sim- 
ulated values of 9 should represent the real distribution 
of the inclination angles in a statistical sense. 

Figure 3 shows the correlation between the simulated 
cos^? and the observed b/a. The whole sample is divided 
into 15 subsets in b/a, with each subset containing a sim- 
ilar number of galaxies. In the figure, points are median 
values while errorbars are the 16-84 percentiles of the 
distribution. The distribution of cos0 is also plotted in 
Figure 1 as the dashed lines. As one can see, cos follows 
a nearly flat distribution as expected. However, it should 



4 



Z. Shao, Q. Xiao, S. Shen, others 



0.8 




b/a 



-20 
M"-51og{h) 



Fig. 3. — The correlation between the apparent axis ratio 6/a and 
the simulated inclination cos 9. Points are located at the median 
value of b/a and cos 9, and error bars represent the 16-84 percentiles 
of the distribution. 



be emphasized that cos 6 is an inclination parameter that 
has only statistical meaning. For individual galaxies, the 
uncertainty in cos6 is as large as the scatter shown in 
Figure 3. 

3. GALAXY LUMINOSITY FUNCTION: DEPENDENCE ON 
INCLINATION 

3.1. Luminosity function estimator 

In estimating the luminosity function (LF) of galax- 
ies, we use the maximum- likelihood method proposed 
by Sandagc, Tammann & Yahil (1979, hereafter STY). 
This method assumes that the LF of galaxies has the 
Schechter form (Schechter 1976), 

(l){L)dL = <j)*{L/L*)" exp{-L/L*)d{L/L*) , (6) 

where L* is the characteristic luminosity, a the faint-end 
slope, and (t>* the overall amplitude of the LF. The STY 
method uses such a luminosity function to calculate the 
probability for a galaxy with luminosity Li and redshift 
Zi to be included in a magnitude- limited sample: 



Pi 



ct>{Li)AL 



(7) 



where imin(-2i) and imax(-2i) are the lowest and highest 

luminosities that a galaxy at redshift Zi can have in or- 
der for it to be included in the sample. The likelihood 
function is defined as 



(8) 



where the product extends over all galaxies in the sample, 
and its maximization provides an estimate for a and L* 
(or equivalently the corresponding absolute magnitude 
M*). 

We estimate the LFs for spiral galaxies with different 
inclinations in all the 5 SDSS bands. In order to make 
meaningful comparisons among the results, only galaxies 



Fig. 4. — The luminosity functions for the nearly face-on sub- 
sample (open circles) and the nearly edge-on subsample (solid cir- 
cles). Smooth curves are the luminosity functions obtained using 
the STY method. Dotted curves show the best fit of the data 
points with a Schechter function. Solid curves are the fit results 
assuming a = (ctsub) = —1.25. Points and errors arc obtained from 
the SWML method. Histograms near the bottom arc the numbers 
of galaxies in bins of absolute magnitude. The results for edge-on 
galaxies arc shifted vertically by a factor of two for easy to compare 
the shape of LFs. 



in the luminosity range from M* + 3.5 to M* — 2.0 are 
used. In practice, an initial value of M* was chosen and 

an iterative procedure was used to determine M* and 
other LF parameters. Some extreme bright galaxies are 
rejected automatically by this selection. 

As comparison, we also estimate each LF using the 
step-wise maximum-likelihood (SWML) method of Efs- 
tathious et al. (1988). Here the LF is represented by 
a non-parametric step function and the maximization of 
the likelihood function is used to determined the relative 
amplitudes at all the steps. 

3.2. Results 

In Table 2 we list the fitting parameters of the LFs 
of spiral galaxies in the 5 SDSS bands, using the STY 

method. SWML method gives similar results; they are 
not presented here in order to save space. The subscript 
in aobs and M*^^ denote that these are derived di- 



'obs' 

rectly from observational data, to distinguish them from 
the quantities for other LFs discussed in this paper. 
Comparing the results with that obtained by Blanton 
et.al. (2001), for all SDSS galaxies, we see that our LFs 
are steeper in the faint end (i.e. a is more negative) in all 
5 bands. This is consistent with the fact that the faint- 
end slope for late-type galaxies is steeper than that for 
early types (Nakamura et al. 2003) . In addition, the typ- 
ical color of galaxies in our sample, as represented by the 
values of M*^^^ in the 5 bands is bluer that for the total 
sample, as is expected from the fact that spiral galaxies 
are on average bluer than early-type galaxies. To study 
the dependence of the LF on galaxy inclinations, we di- 
vide sample of each band into 15 subsamples according 
to the apparent axis ratio b/a in the r band (see § 2.2 



TABLE 2 

Parameters of Luminosity Function pgr Spiral Galaxies in 
Each of the 5 SDSS Bands 



Band 


-^spiral 


CKobs 




u 


23757 


-1.47±0.02 


-18.32±0.02 


9 


42620 


-1.36±0.01 


-19.60±0.01 


r 


60523 


-1.34±0.01 


-20.26±0.01 


i 


51593 


-1.36±0.01 


-20.62±0.01 


z 


43030 


-1.31±0.01 


-20.74±0.01 



and Figure 3 for details). LF is estimated for each of the 
subsamples in each of the 5 SDSS wave bands. As an 
example, Figure 4 shows the LFs for two extreme sub- 
samples, one for nearly cdgc-on galaxies, and the other 
for nearly fkcx}-on galaxies, in the r-band. Note that the 
results obtained from the SWML method (open and solid 
circles) match extremely well with those obtained from 
the STY method (solid lines) , showing that the Schechter 
function is a valid assumption for the LFs for spiral galax- 
ies with different inclinations. Clearly, both subsamples 
have similar a, but edge-on galaxies are systematically 
fainter than face-on galaxies, by about ~ 0.65 magni- 
tude. In Table 3, we list the LF parameters obtained 
from the SYT method for all the inclination subsamples, 
and we use the subscript 'sub' to denote the parame- 
ters for subsamples. The changes of asub and M*^y^ with 
b/a are also plotted in Figure 5 and in the left panel 
of Figure 6, respectively. There is a systematic change 
of M*^jj, with face-on galaxies having a brighter charac- 
teristic magnitude than edge-on galaxies. This change 
of M*^^ with inclination becomes smoother if we fix the 
faint-end slope a to be the average value of the 15 sub- 
samples (dotted lines in Figure 5), as one can see from 
the right panel of Figure 6. These values are denoted 
as M*^^ 2 and also listed in Table 3. In addition, the 
inclination-dependence of M*^^ is stronger in blue band, 
as is expected if the dependence is caused by dust ex- 
tinction. 

A systematic change of a with 6/a is expected if dust 
extinction depends on galaxy luminosity for L < L* . 
The absence of any strong systematic dependence of a 
on b/a in the data suggests that dust extinction does 
not depend on galaxy himinosity strongly, at least for 
sub-L* galaxies. Note, however, that the dust extinction 
considered here is relative to that for face-on galaxies. 
If dust extinction in face-on galaxies depends strongly 
on luminosity, then a strong luminosity-dependent dust 
extinction cannot be excluded. We will discuss this in 
more detail in S 5.1. 
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Fig. 5. — The dependence of the faint-end slope of the LF, a, on 
the axis ratios of galaxies {b/a). Dotted lines for (osub) of table 3; 
Solid lines for a^hs of table 2; Dashed line for ao of table 5 (see 
§ 4.1). 
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Fig. 6. — The dependence of the characteristic magnitude of the 
LF, M*, on the axis ratios of galaxies {b/a). In the left panel, 
the values of M* arc taken from the fit with both M* and a as 
free parameters. In the right panel, the values of M* are obtained 
from the fit assuming a = (ogub)- Solid lines for M*^^^ of table 2; 
Dashed line for Mo of table 5 (see § 4.1). 
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TABLE 3 

Luminosity Function Parameters for Spiral Galaxies with 
Different Inclinations. 



xJCJii-i-Vj. J. y SUD, 1 


yyj 1 J mccl 


( COS 




su 


K 1 

-'^-^sub.l 


su ,2 


(asub) 




U 1256 


0.21 


0.09 


-1 


21±0.09 


-17.00±0.06 


1260 


-1 


27 


-17.04±0.06 


1216 


0.28 


0.17 


-1 


14±0.08 


-17.27±0.06 


1223 






-17.36±0.07 


1373 


0.34 


0.24 


-1 


23±0.08 


-17.65±0.06 


1374 






-17.68±0.06 


1440 


0.39 


0.29 




09i0.07 


-17.65i0.06 


1451 






-17.79±0.06 


1519 


0.44 


0.36 




32±0.07 


-17.99±0.07 


1516 






-17.95±0.06 


1544 


0.49 


0.42 


-1 


23±0.07 


-17.99±0.06 


1545 






- 18.02±0.06 


lb58 


0.55 


0.49 




46zb0.07 


- 18 . 27±0 . 07 


1657 






- 18. 15±0.06 


1604 


0.59 


0.54 




34±0.07 


-18.28±0.06 


1606 






-18.24±0.07 


1672 


0.64 


0.60 




33i0.07 


-18.33i0.06 


1672 






-18.30i0.06 


1685 


0.69 


0.66 




33±0.07 


-18.35±0.06 


1685 






-18.31±0.05 


1665 


0.74 


0.72 




47±0.07 


- 18.56i0.07 


1664 






-18.41i0.07 


1750 


0.79 


0.77 




32±0.07 


-18.42±0.06 


1746 






-18.38±0.05 


1709 


0.83 


0.82 


-1 


34±0.07 


-18.41±0.06 


1705 






-18.36±0.06 


1765 


0.88 


0.88 


-1 


15 ±0.06 


-18.44±0.06 


1769 






- 18.52±0.06 


1860 


0.94 


0.94 




10i0.06 


-18.42±0.05 


1870 






-18.54±0.05 


Q 2672 


0.21 


0.09 


— 

-1 


25±0.05 


-18.65±0.05 


2671 


-1 


24 


-18.64±0.05 


2598 


0.28 


0.17 


-1 


25±0.05 


-18.87±0.05 


2598 






-18.86±0.05 


2726 


0.34 


0.24 


-1 


30±0.05 


-19.27±0.05 


2726 






-19.23±0.05 


2788 


0.39 


0.29 




24zb . 05 


-19. 22 ± . 04 


2788 






- 19.22±0.05 


2841 


0.44 


0.36 




33i0.05 


-19.44±0.05 


2842 






-19.38±0.04 


2834 


0.49 


0.42 


-1 


29±0.05 


-19.44±0.05 


2834 






-19.41±0.05 


2894 


0.55 


0.49 




25±0.05 


-19.53±0.05 


2897 






-19.52±0.04 


2853 


0.59 


0.54 




25±0.05 


-19.56±0.05 


2852 






-19.56±0.05 


2887 


0.64 


0.60 




26±0.05 


-19.65±0.05 


2889 






-19.63±0.04 


2921 


0.69 


0.66 




32zh0.05 


-19.67±0.05 


2920 






-19.61±0.05 


2897 


0.74 


0.72 


"""" 


31zb0.05 


-19. 74±0.05 


2896 






-19.69±0.04 


2925 


0.79 


0.77 




23i0.05 


- 19.74i0.05 


2926 






-19. 75i0.05 


2839 


0.83 


0.82 




22±0.05 


-19.66±0.05 


2838 






-19.68±0.04 


3008 


0.88 


0.88 




09i0.05 


- 19.63i0.04 


3014 






-19. 74i0.04 


2981 


0.94 


0.94 




02±0.05 


-19.66±0.04 


2983 






-19.82±0.04 


r 4093 


0.21 


0.09 


— 

-1 


25±0.04 


-19.50±0.04 


4093 


-1 


25 


-19.50±0.04 


4075 


0.28 


0.17 


-1 


27±0.04 


-19.68±0.04 


4075 






-19.67±0.04 


4047 


0.34 


0.24 


-1 


26±0.04 


-19.99±0.04 


4048 






-19.98±0.04 


4026 


0.39 


0.29 




25i0.04 


-20.00i0.04 


4027 






-20.00i0.04 


4011 


0.44 


0.36 




32i0.04 


-20.20i0.04 


4010 






-20. 16i0.04 


4092 


0.49 


0.42 




31±0.04 


-20. 16i0.04 


4089 






-20. 1 li0.04 


3987 


0.55 


0.49 


"""" 


30zb 0. 04 


-20. 3 1 ± 0. 04 


3988 






-20.28 ±0.04 


3987 


0.59 


0.54 


"""" 


27±0.04 


-20.28±0.04 


3987 






-20. 27±0. 04 


4072 


0.64 


0.60 


-1 


29±0.04 


-20.34±0.04 


4071 






-20.32±0.04 


4042 


0.69 


0.66 




36±0.04 


-20.40±0.04 


4043 






-20.32±0.04 


4030 


0.74 


0.72 




32±0.04 


-20.39±0.04 


4033 






-20.34±0.04 


4061 


0.79 


0.77 




18±0.04 


-20.32i0.04 


4059 






-20.37±0.04 


3979 


0.83 


0.82 




21±0.04 


-20.29i0.04 


3983 






-20.33i0.04 


4081 


0.88 


0.88 




17zt0.04 


-20.34±0.04 


4085 






-20.40±0. 04 


^ 4014 


0^94 


0^94 




07zb 0. 04 


-20.31 ±0.04 


4016 






-20.44±0.04 


i 3700 


0.21 


0.09 


— 

-1 


29±0.04 


-20.00±0.04 


3700 


-1 


30 


-20.00±0.04 


3603 


0.28 


0.17 


-1 


23±0.05 


-20.08±0.04 


3604 






-20.13±0.04 


3585 


0.34 


0.24 


-1 


30±0.05 


-20.43±0.04 


3585 






-20.43±0.05 


3543 


0.39 


0.29 




40±0.05 


-20.52±0.05 


3546 






-20.43±0.04 


3456 


0.44 


0.36 




31±0.05 


-20.57±0.05 


3456 






-20.56±0.04 


3503 


0.49 


0.42 




44zh0.05 


-20.62±0.05 


3500 






-20.52±0.05 


3394 


0.55 


0.49 


"""" 


34zb 0. 05 


-20. 71±0.04 


3393 






-20.68±0.05 


3364 


0.59 


0.54 




30i0.05 


-20.67±0.05 


3365 






-20.67±0.05 


3401 


0.64 


0.60 




27i0.05 


-20.64i0.04 


3402 






-20.67i0.04 


3332 


0.69 


0.66 




39±0.05 


-20.77±0.05 


3336 






-20.70±0.05 


3355 


0.74 


0. 72 


-1 


37zb0.05 


-20. 74±0.05 


3355 






-20.69±0.05 


3385 


0.79 


0. 77 


-1 


17zb0.05 


-20.61±0.04 


3387 






-20.70±0.04 


3270 


0.83 


0.82 


-1 


19±0.05 


-20.58±0.04 


3269 






-20.65±0.04 


3381 


0.88 


0.88 




28±0.05 


-20.70±0.05 


3379 






-20.71±0.05 


3338 


0^94 


0^94 




18i0.05 


-20.68i0.04 


3335 






-20. 77i0.04 


Z 3275 


0.21 


0.09 


— 
-1 


32±0.05 


-20.35±0.05 


3273 


-1 


27 


-20.31±0.05 


3229 


0.28 


0.17 


-1 


21±0.05 


-20.31±0.05 


3233 






-20.34±0.04 


3147 


0.34 






22zb0. 05 


-20. 59±0. 05 








-20. 63 ±0.05 


3086 


0.39 


0.29 


-1 


42 ±0.05 


-20.72±0.05 


3076 






-20, 59 ±0.05 


2928 


0.44 


0.36 


-1 


22±0.06 


-20.70±0.05 


2929 






-20.74±0.05 


2934 


0.49 


0.42 


-1 


39±0.06 


-20.77±0.05 


2932 






-20.69±0.05 


2859 


0.55 


0.49 


-1 


35±0.06 


-20.86±0.05 


2857 






-20.80±0.05 


2724 


0.59 


0.54 


-1 


22±0.06 


-20.77±0.05 


2724 






-20.80±0.05 


2775 


0.64 


0.60 


-1 


20±0.06 


-20.76±0.05 


2778 






-20.81±0.05 


2755 


0.69 


0.66 


-1 


40±0.06 


-20.87±0.05 


2752 






-20.78±0.05 


2640 


0.74 


0.72 


-1 


23±0.06 


-20.77±0.05 


2641 






-20.80±0.05 


2704 


0.79 


0.77 


-1 


16±0.06 


-20.75±0.05 


2706 






-20.81±0.05 


2603 


0.83 


0.82 


-1 


24±0.06 


-20.73±0.05 


2604 






-20.75±0.05 


2682 


0.88 


0.88 


-1 


28±0.06 


-20.81±0.05 


2684 






-20.81±0.05 


2705 


0.94 


0.94 


-1 


17±0.06 


-20.79±0.05 


2706 






-20.85±0.05 



Samples of all 5 bands arc divided into 15 sub-samples by b/a, 
while (6/a)nicd and (cos(?)jncd show the median values of inclina- 
tion of each subsamplc. (osubi ^s*ub l) denoted the fitting param- 
eter of LF for each sub-sample, and M*^^ ^ are fitting values of 
characteristic magnitude when we fixed the slope parameter to be 
the average value of 15 sub-samples, (agub). 
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Flc;. 7. — The inclination-dependence of M* (M*^^ ,^ listed in 
table 3, see text for details) and its best linear fit. Solid and open 
circles are used to distinguish two successive wave-bands. 



4. INTERNAL DUST EXTINCTION 

Having shown that the LF of spiral galaxies depends 
systematically on the inclination, wc now model such de- 
pendence in terms of dust extinction. 

4.1. Empirical extinction models 

As a simple model for the dust extinction, wc assume 
that the change in magnitude of a galaxy. AM, due to 
dust extinction is proportional to log(6/a), so that 

AMi{b/a) = M{b/a) - Afi(l) = -71 log(6/a) , (9) 

where 71 parameterizes the amplitude of the dust extinc- 
tion relative to the face-on value (e.g. Giovanelli et al. 
1994; Masters, Giovanelli & Haynes 2003). As discussed 
in § 2.2, a quantity that may better describe the inclina- 
tion is the corrected cosine of the inclination angle, cos^. 
We therefore consider another model in which 

AM2(cos6') = M(cos6') -M2(l) = -72 log(cos 6*) . (10) 

In the above expressions, the subscripts, '1' and '2' de- 
note quantities in the two models. Note that Afi(l) and 
M2(l) arc the magnitudes corresponding to b/a = 1 and 
cos6' = 1, i.e. face-on disk. 

Since the inclination-dependence of asub is not strong, 
the main effect of dust extinction is to change M*^^^. The 
change of M*^^^^ with inclination may be used to infer 
the average dust extinction for all spiral galaxies in the 
sample. The results of the least square fits of equations 
(9) and (10) to the observed M*^^ .2-log{b/a) and M*^^ .^- 
log(cos 9) relations are shown in Fig. 7, while the val- 
ues of the best fit 7 and Af*(l) arc listed in Table 4 
along with the corresponding values. There is signifi- 
cant difference between the results based on the apparent 



TABLE 4 

Fitting parameters of the linear relationship between M* 

AND THE LOGARITHMIC OF THE INCLINATION. AlSO SEE SOLID 

LINES IN Figure 7. 



Baud 




Ti 






~|2 


xi 


u 


-18.65±0.03 


2.19±0.08 


18.82 


-18.58±0.02 


1.48±0.06 


7.24 


9 


-19.90±0.02 


1.68±0.06 


31.36 


-19.84±0.02 


1.14±0.04 


15.31 


r 


-20.53±0.02 


1.37±0.0.5 


48.28 


-20..18±0.02 


0.92±0.04 


26.29 


i 


-20.84±0.02 


1.08±0.06 


41.61 


-20.81±0.02 


0.73±0.04 


26.77 


z 


-20.92±0.02 


0.80±0.07 30.98 


-20.89±0.02 


0.54±0.04 


23.71 



axis ratio hja and those based on the corrected inclina- 
tion angles Q. For all of the 5 wave bands, xh < xfi 
suggesting that a linear model works better for the cor- 
rected inclination angle. The values of 72 are all smaller 
than 71, because the relation between bja and cos^ has 
a slope that is different from 1 (see Figure 3). However, 
the ratios between the 7 values in different wave bands 
are quite similar in the two models. As we will see in 
§ 4.2, it is these ratios that describe the dust extinction 
curve. In the maximum-likelihood estimate of the LF 
used here, dust extinction can be incorporated in a more 
elegant way. Since equation (10) provides a reasonable 
description of the mean relation between the dust extinc- 
tion and the inclination angle a similar relation may 
be used for individual galaxies. In this case, the luminos- 
ity function of spiral galaxies with observed luminosity 
L and inclination angle Q can be parameterized by 

0(L,0)oc [(L/L*)(cos0)-°-^'^]"exp[-(L/L*)(cos0)-°-4^] , 

(11) 

Thus, the value of 7 (assumed to be independent of 
galaxy luminosity) can also be determined through the 
maximum likelihood method (either the STY method or 
the SWML method) that determines the LF. We have 
carried out such a maximum likelihood analysis using 
the STY method. Again, an iterative procedure is used 
so that only galaxies with luminosities in the range be- 
tween M* -h 3.5 and M* - 2.0 are used in the fitting. 
Since cos^ for individual galaxies are obtained statisti- 
cally from Monte-Carlo simulations, different realizations 
may give different results. We therefore construct 25 dif- 
ferent realizations and obtain the values of a, M* and 
7 for each of them. The average values of a, M* and 
7 obtained in this way are listed in Table 5 as ao, 
and 73, where the subscript '0' is used to denote 'face-on' 
galaxies, and the subscript '3' in 73 is used to distinguish 
it from 71 and 72 obtained from the M^^^ .^r^o^i^ja) rela- 
tion and the M*^^^ 2"log(cos ff) relation. The uncertainties 
on these parameters include both fitting errors and the 
scatter among the 25 realizations. The error bars are all 
very small, because each sample is sufficiently large to 
represent the distribution of cos^ faithfully. Note that 
the values of 73 are similar to that of 72, indicating that 
the method based on M* and that on individual galaxies 
give the same results. 

4.2. Extinction curves 

The wavelength-dependence of 7 can be used to con- 
strain the extinction curve for spiral galaxies. As a simple 
model, we assume that the effective extinction curve in 
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TABLE 5 

Corrected Luminosity Function (eq. 11) for Spiral 
Galaxies in Each of the SDSS Band 



Band N^pira.! 


ao 


M* 


73 


u 


23380 


-1.28±0.02 


-18.60±0.02 


1.59±0.05 


9 


42150 


-1.22±0.01 


-19.89±0.01 


1.24±0.03 


r 


60100 


-1.25±0.01 


-20.54±0.01 


0.92±0.02 


i 


51290 


-1.28±0.01 


-20.87±0.01 


0.80±0.02 


z 


42820 


-1.26±0.01 


-20.96±0.01 


0.65±0.03 



. These values represent the parameters for pure face-on spiral 
galaxies. 



TABLE 6 

Fitting Results for the Extinction Curve (eq. 13). 



Case 1 (using 7i) 1.45±0.02 0.97±0.07 3.94 
Case 2 (using 72) 0.98±0.02 0.97±0.07 3.71 
Case 3 (using 73) l.OSiliO.Ol 0.96±0.04 1.50 



the optical has a power-law form, 

rA(cos6l) = ry(cos6')(A/5500l)"" , (12) 

where A is the wavelength, t\ is the optical depth at 
A, Ty is the optical depth at the V^-band (centered at 
A = 55OOA) and is a function of inclination angle, and 
n is the power index describing the shape of the dust 
extinction curve. If wc neglect the difference in the dust 
extinction curve for different galaxies and compare the 
dust extinction at given inclination angle but in different 
wavelength, then 7 is proportional to r. Thus, we can 
write 

7A = 7y(A/5500A)-", (13) 

where jv is the value of 7 at 5500A. 

We use a least-square fitting to estimate "/v and n from 
the values of 7 in the 5 SDSS bands, and the results are 
listed in Table 6 for the three cases. Cases 1, 2 and 3, 
corresponding to the use of 71, 72 and 73, respectively. 
In Figure 8, wc plot the fitting results along with the 
data of 7 used in the fit. As one can sec, the power- 
law model is a good assumption in all cases. It is also 
interesting to note that, although the values of 71, 72 
and 73 are different, the power-law indices obtained from 
them arc very similar, with n = 0.96 ± 0.04. (12) we 
obtained from the SDSS data are extrapolated to the 
near-infrared, and compared with the data obtained by 
Masters et al. (2003) for spiral galaxies with log(a/6) < 
0.5 in the 2MASS: 77 = 0.48 ± 0.15, jh = 0.39 ± 0.15 
and jKs = 0-26 ± 0.15. In addition, we also plot results 
obtained by others in various bands. TuUy et al. (1998) 
gave 7b = 1.0 -1.5 for different types of spiral galaxies. 
Courteau (1996) obtained Jrcunn = 0-95 using 349 Sb-Sc 
UGC galaxies with 0.27 < log{a/b) < 0.70. Han (1992) 
found 77 to be in the range from 0.51 to 0.91 for different 
types of galaxies in the / band. Finally, GiovancUi ct al. 
(1994) obtained 7/ = 0.95 - 1.15. Remarkably, all of the 
results follow well the power-law we obtained from the 
SDSS data. 




o 7i (Giovanelli et al. 1994) 

« 7i (Han 1992) 

^ "/j.H.Ks (Masters et al. 2003) 



" 5000 10" 

A(A) 

Fic;. 8. — The best fit of tfic extinction curve in optical wave- 
bands. The values of 7 are those given by the different methods 
described in § 4.2. Data points show the values of 7 obtained by 
other investigators. 



It should be pointed out, however, that the dust ex- 
tinction we obtained here is the effective extinction for an 
entire galaxy, rather than that from some specific regions 
within a galaxy. Such extinction depends not only on the 
properties of the dust grains but also on the distribution 
of dust in individual galaxies. The index, n = 0.96±0.04, 
obtained is shallower than that obtained by observations 
of the Milky Way, the LMC and the SMC, n = 1.1 - 1.5, 
but steeper than the value n = 0.7 obtained by Char- 
lot et al. (2000) based on the assumption that dust in 
a spiral galaxy has a patchy distribution. Our result of 
n ~ 1 implies that dust distribution in spiral galaxies on 
average are not as patchy as assumed in Chariot et al. 
(2000). 

5. EXTINCTION-CORRECTED LUMINOSITY FUNCTION 

5.1. Luminosity- dependence of extinction 

In the analysis presented above, we have assumed that 
the value of 7 is independent of luminosity. In this sub- 
section we examine whether or not this assumption is 
valid. To do this, we consider a model in which 7 is a 
linear function of the magnitude, 

j{M) = j*[l + P{M-M*)]. (14) 

We treat (3 as the fourth free parameter in fitting the 
LF. The result shows that has a very small value so 
that the difference between 7 and 7* is less than 5%. 
Furthermore, if a is kept constant in the fitting, f3 is 
even smaller. Thus, the luminosity-dependence of 7 is 
weak and will be neglected in the following. Note, how- 
ever, that 7 in our definition only characterizes the ex- 
tinction at a given inclination relative to that of face-on 
galaxies. Thus, even if 7 is independent of galaxy lu- 
minosity, the absolute value of the extinction may be 
luminosity-dependent, if the extinction of face-on galax- 
ies depends on luminosity. Unfortunately our method, 
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which is based on a comparison between inclined galax- 
ies with face-on galaxies, cannot be used to obtained the 
absolute value of the extinction. Because of this, here 
we use an independent set of observational data to probe 
the luminosity-dependence of extinction and to exam- 
ine how such data can be used together with our results 
to obtain an extinction-corrected luminosity function for 
spiral galaxies. 

Kauffmann et al. (2003) estimated the dust extinc- 
tion in the SDSS z-band, Az, for galaxies in the SDSS 
DR2, using the 4000 A break strength and the Balmer ab- 
sorption line index HSa to constrain the star formation 
history and then estimating the dust extinction from the 
difference between the model and observed g — r and r — i 
colors. We combine their data with our sample. It should 
be noted that the 4GGGA breaks and the USa indices are 
estimated from the SDSS fiber spectra within a 3-arcsec 
aperture. Thus, the values of given by Kauffmann et 
al. may be biased towards the central parts of galaxies, 
especially for nearby galaxies with large apparent sizes. 
Figure 9 shows A^ versus luminosity for galaxies in our 
sample. The small dots in the upper panel show all galax- 
ies in the 'face-on' subsample with 0.86 < b/a < 0.91 . 
Open circles with error bars represent the median and 
their errors within given bins of M^- As a comparison, 
the solid circles show the median values of Az for 'edge- 
on' galaxies with 0.25 < b/a < 0.31. The solid circles 
are shifted by 0.5 magnitude towards the bright end ac- 
cording to the difference in M*^^^ for these two subsam- 
ples. It is clear that dust extinction does depend on 
luminosity for both face-on and edge-on galaxies. The 
dependence is not a simple nionotonic; relation. The ex- 
tinction increases with luminosity at the faint part and 
decreases slightly towards the bright end. The behav- 
ior at the faint end may indicate that the dust opacity 
is larger in brighter galaxies. The slight decrease of Az 
with luminosity at the bright end may be due to the fact 
that the contribution of the central bulge, which may be 
less dusty, becomes more important for iDrighter galaxies. 
Note that the mean Az-Mz relation for the face-on sub- 
sample is roughly parallel to that for the edge-on subsam- 
ple, except perhaps for faint galaxies (see the lower panel 
of Figure 9). Since "fz oc AA^, this result implies that 7^ 
is quite independent of luminosity, consistent with our 
result based on the luminosity functions. 

The difference in Az between face-on and edge-on 
galaxies is about 0.4 magnitude, in good agreement with 
the difference in M*^^^ (0.5 mag) between these two sub- 
samples. We have also checked the M^-dependence of Az 
for other subsamples of b/a, and found that the relations 
are all roughly parallel to each other, with the relative 
amplitudes very similar to those obtained from fitting 
the LFs (see Table 3). This agreement is remarkable, 
because the A^ values given by Kauffmann et al. (2003) 
are based on totally different considerations. 

5.2. Extinction- corrected luminosity function 
5.2.1. Luminosity function of face-on galaxies 

Based on the results presented above, we can have sev- 
eral estimates for the LF of face-on spiral galaxies. The 
first one is based on (asuh, M*^^ -^^), or, ((a'sub)- A^*ib,2): 
obtained from the subsamples with big values of axis ra- 
tio (see lines 14 and 15 in Table 3). The second is given 
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Fig. 9. — The luminosity-dependence of the dust extinction in 
the x-band (Az). In the upper panel, small dots axe the values of 
Az obtained by Kauffmann et. al. (2003) for individual fa<;e-on 
spiral galaxies with 0.86 < b/a < 0.91. Open circles with error 
bars represent the median value of Az for a given Mz ■ Solid circles 
arc the same as open ones but for galaxies with 0.25 < b/a < 
0.31 (almost edge-on). The lower panel shows the difEcrciice of 
Az between face-on and edge-on galaxies. Solid lines here show a 
fitting by eye of the data points. 



by ((asub),Mr(l)) or ((a,ub), M^l)), where Mi*(l) and 
M|(l) are the values of M* for face-on galaxies obtained 
from fitting the M*j^(,-inclination relations (see Table 4). 
The most rigorous estimate may be given by the values 
of (aoi-^o) Table 5, obtained through the maximum 
likelihood analysis that includes 7 as the 3rd parameter 
(eq. 11). Note that for all the 5 SDSS bauds, the values of 
(Qsub) and ao are similar, but both are slightly less neg- 
ative than aobs, the faint end slope of directly measured 
LF of spirals (§ 3.2). Furthermore, unlike ctobs, which 
has more negative values in the bluer bands, both (agub) 
and ao are quite independent of wave-bands. This may 
be explained by the fact that dust extinction not only 
reduces L*, but also steepens the LF for spiral galaxies. 
We have tested this effect using a Monte Carlo simula- 
tion. We generated a sample of galaxies with a given LF, 
and with each galaxy assigned a random orientation. We 
then made each galaxy dimmer according to its inclina- 
tion angle. We found that the resulting LF is steeper 
than the original LF. 

The corrections in the characteristic magnitudes are 
also significant. For face-on galaxies, Mj*(l) , M|(l) and 
Mq are all brighter than M*y^^ for about 0.2-0.3 magni- 
tudes. Figure 10 shows a comparison of M*s for all the 
5 bands. According to the extinction curve (t oc A""'^), 
the change in M* is expected to be larger for bluer band. 
That such dependence is not secin clearly in Figure 10 is 
largely caused by the change in a and the degeneracy 
between a and M*. 
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Fig. 10. — A comparison between the values of M* obtained from 
different cases in different bands. LFobg labels the results of the 
observational sample without any dust correction. LFq labels the 
results for face-on galaxies obtained by correcting the inclination 
effects. LFint labels the results for 'extinction-corrected' luminosity 
functions. 



5.2.2. Intrinsic Luminosity Function 

According to Figure 9, the values of for face-on 
galaxies are typically 0.3-0.4 magnitudes. In principle, 
we can use the values of for individual galaxies and 
estimate a dust-corrected LF. The LFs in other bands 
can also be obtained using the the dust extinction curve. 
Unfortunately, this approach is impractical, because the 
uncertainty in A^ is quite large {AA^ > 0.2) and can sig- 
nificantly broaden the LF. On the other hand, for galax- 
ies with a given luminosity, the true values of Az may 
have intrinsic scatter, which should be taken into ac- 
count in the dust correction. Because of this uncertainty, 
here we consider two extreme cases. In one, we use the 
Az values given by Kauffmann et al. (2003) to make 
correction for each galaxy, and estimate the LF for the 
'corrected' sample. The LF parameters obtained in this 
way are a = -1.43 ± 0.01 and M* = -21.52 ± 0.01. In 
the second case, we average the values of A^ for galaxies 
with similar inclination angles, and use them to correct 
individual galaxies. The LF parameters obtained in this 
way are a = -1.28 ± 0.01 and M* = -21.27 ± 0.01. As 
discussed above, the true LF is expected to be between 
these two. Our Monte Carlo simulation showed that the 
uncertainty in A^, about 0.2, can account for most of 
the difference between these two results. So the second 
estimation of LF may close to the real case. 

For the other 4 bands, direct estimates of the dust 
extinction are not available. We therefore use an as- 
sumption of the dust extinction curve together with the 
average values of A^ (as a function of b/a) to make dust 
corrections for individual galaxies. We list the param- 
eters for the LFs obtained for the corrected samples in 
Table 7 as aint and Mj*^. For comparison, results are 
given for two extinction curves, one with n = 1.0, as is 
obtained in this paper, and the other is n = 0.7, as given 
in Chariot et. al. (2000). 

In Figure 10, we plot M* in all the five bands for dif- 
ferent cases: the directly observed LF, the LF for face- 
on galaxies, and the intrinsic LFs obtained above. It 
is clear that the differences between different cases are 
quite large. In particular, M* for the corrected LF could 
be 0.5-1.2 magnitudes brighter than the uncorrected one 



TABLE 7 

Parameters of Luminosity Function for Spiral Galaxies 
with correction of dust extinction. 



n=1.0 


n=0.7 


Band 


-^spiral 


Oint 




^ spiral 
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22887 


-1.46±0.01 


-19.58±0.02 


22902 


-1.50±0.01 


-19.28±0.02 
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41587 


-1.32±0.01 


-20.56±0.01 


41582 


-1.32±0.01 


-20.32±0.01 


r 


59545 


-1.26±0.01 


-20.99±0.01 


59536 


-1.26±0.01 


-20.90±0.01 


i 


50718 


-I.30it0.01 


-21.23±0.01 


50719 


-1.30±0.01 


-21.19±0.01 


z 


42311 


-1.28±0.01 


-21.27±0.01 


42311 


-1.28±0.01 


-21.27±0.01 



from 2 to u bands. This suggests that, without correcting 

internal extinction of galaxies, the L* for spiral galaxies 
can be seriously undcr-cstimatcd. 

6. SUMMARY AND DISCUSSION 

Using a samples of 61506 spiral galaxies selected from 
the SDSS DR2, we study the luminosity functions of spi- 
ral galaxies with different inclination angles. The ap- 
parent axis ratio, 6/a, is used as an observational incli- 
nation indicator to define subsamples of spiral galaxies 
at different inclinations, and we use a Monte Carlo pro- 
cess to connect hj a to the corrected inclination angle, 0, 
disk thickness, and ellipticity e. There is a systematic 
change of the LF with inclination angle: the character- 
istic luminosity L* decreases with increasing inclination 
angle, while the faint-end slope depends only weakly on 
inclination. 

The inclination-angle dependence of the galaxy lumi- 
nosity function is consistent with the expectation of a 
simple model where the optical depth is proportional to 
cos^, and we use a likelihood method to recover both 
the amplitude of the extinction (relative to the face-on 
value), 7, and the luminosity function for galaxies viewed 
face-on. We found that the results obtained from dif- 
ferent methods are all consistent with each other, and 
the characteristic magnitude for face-on spirals is about 
0.2 ~ 0.3 magnitudes brighter than the average popula- 
tion. 

We found that the value of 7 is quite independent of lu- 
minosity in a given band. The values of 7 obtained in this 
way for the 5 SDSS bands are used to constrain the shape 
of the extinction curve, assuming t\ = ry (A/5500j4.)~". 
We find n ^ 0.96 ± 0.04. 

Using the z-band dust extinction given by Kauffmann 
et al. (2003), together with the inclination-dependence of 
the LF we obtained, we derive an 'extinction-corrected' 
luminosity function for spiral galaxies. Dust extinction 
makes a significant change in M*, and the characteris- 
tic luminosity of the 'dust-corrected' LF is about 0.5 to 
1.2 magnitude brighter than the uncorrected LF from 
the z- to the u-bands. This suggests that the luminosity 
function of spiral galaxies may be significantly underes- 
timated in blue bands, if internal dust correction is not 
made. This may have important implications when com- 
paring model predictions of the luminosity function with 
observations. 

As mentioned in the introduction, the dimming in lu- 
minosity from facc-on to edge on is expected if galaxy 
disks arc optically thick. Our results therefore give 
support to the assumption that most of the disks in 
our galaxy sample are optically thick. Note that the 
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inclination-dependent dimming exists not only for faint 
galaxies but for galaxies over the entire luminosity range. 
Note also that we are using a well-defined flux-limited 
sample and the selection bias is taken into account in 
the LF estimate. It is therefore unlikely that the effect 
we find here is due to the magnitude limit in the sam- 
ple. Systematic effects of surface brightness relative to 
the inclination is also unlikely to play an important role 
here. As discussed in Blanton et al. (2001), the surface- 
brightness limit used in the SDSS has negligible effects on 
the LF at Mr < —18, and so the inclination-dependence 
of the LF at Mr < —18 cannot be explained as due to 
the miss of low-surface face-on galaxies. We therefore 
conclude that our results are best explained by assuming 
that bright disks are optically thick. 
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